f-Sum Rule and Unconventional Spectral Weight Transfer in Graphene 
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We derive and analyze the f-sum rule for a two-dimensional (2D) system of interacting electrons 
whose behavior is described by the Dirac equation. We apply the sum rule to analyze the spectral 
weight transfer in graphene within different approximations discussed in the literature. We find 
that the sum rule is generically dominated by inter-band transitions while other excitations produce 
sub-leading behavior. The f-sum rule provides strong constraints for theories of interacting electrons 
in graphene. 
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I. INTRODUCTION. 

Exact sum rules provide a very insightful way to test 
the validity of certain approximations in the theory of the 
interacting electron liquid*^ In those systems, electrons 
interact via Coulomb forces, which are long-ranged, mak- 
ing perturbation theory questionable. A common way to 
deal with this issue is to perform an infinite diagrammatic 
resummation, such as the random phase approximation 
(RPA), that leads to the screening of the long range part 
of the interaction, effectively reducing it. The result can 
be used to obtain electron response functions, where the 
different excitations arising in the effective theory show 
particular features. This procedure is at the heart of the 
Landau Fermi liquid theory of He 3 and ordinary metals. 2 
In this context, the so-called f-sum rule provides a way 
to check if a particular resummation or approximation 
is describing the relevant excitations, and their relative 
importance.- 

The f-sum rule has been widely used in the context 
of 3D and 2D Landau-Fermi liquids, where interactions 
are switched on adiabatically from an initial gas of free 
Schrodinger electrons. In those systems, in the high den- 
sity regime, the relevant excitations are correctly cap- 
tured by the RPA resummation, leading to Landau quasi- 
particles and a collective plasmon mode. In the RPA ap- 
proximation, perturbation theory can be cured once it is 
realized that the Coulomb potential is screened by the 
electrons. 

Recently, graphene has attracted a lot of attention due 
to the fact that its lattice geometry gives rise to Dirac- 
like electrons described by the 2D Dirac equation 3 . This 
purely 2D crystal was isolated and characterized in the 
laboratory in 2005^ In contrast to an ordinary Fermi Liq- 
uid, neutral or undoped graphene is a semimetal, whose 
density of states vanishes at the Dirac point. In this way, 
the electrons have very poor screening properties, and 
the picture of quasiparticles composed of weakly interact- 
ing electrons could in principle fail. On the other hand, 
experimentally it is known that undoped graphene carri- 
ers behave as non-interacting Dirac electrons, so the role 



played by interactions in graphene is still unclear. Per- 
turbative renormalization group (RG) arguments, which 
are equivalent to the RPA approximation, indicate that 
unscreened Coulomb interactions are marginally irrele- 
vant, only leading to a logarithmic renormalization of 
the Fermi velocity^ Nevertheless, this approach has been 
criticized recently from different directions. In Ref. @, it 
has been argued that new collective modes, arising from 
non-trivial vertex corrections beyond RPA, may play an 
important role in the low-energy description of graphene. 
A different result is obtained when interactions are stud- 
ied in the tight-binding description of graphene. By us- 
ing a ring diagram approximation,— it has been showed 
that electrons in undoped graphene would behave essen- 
tially like a Fermi liquid, in contrast to the marginal 
Fermi liquid behaviour that arises from the RG calcu- 
lation. Therefore, being yet a controversial issue, it is 
important to provide new tools to study the interacting 
problem in graphene, so as to reanalyze the relevance of 
the different electronic excitations in this system. 

In this paper we study the different approximations to 
electron interactions in doped and undoped graphene by 
means of the f-sum rule. In section II we review those 
approaches, putting special emphasis in the different ex- 
citations that are present. In section III we introduce 
the f-sum rule for 2D Dirac electrons and the constraints 
that imposes on approximations to the interacting prob- 
lem. In section IV we study the spectral weight transfer 
in terms of the interaction strength. Section V contains 
the analysis of the results, while in section VI we summa- 
rize the main highlights of the work. We have included 
appendix [A] which contains a detailed derivation of the 
f-sum rule in different situations. 



II. THE THEORY OF INTERACTING 
ELECTRONS IN GRAPHENE. 

Due to the symmetry properties of the honeycomb lat- 
tice there are two inequivalent points in the Brillouin 
zone, K and K'. In the low-energy limit an expansion 
around those points gives an effective Dirac equation rul- 
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ing the motion of the electrons. In dealing with long 
range interactions, the scattering is dominated by small 
momentum transfer allowing us to concentrate on each 
individual point and neglect inter-valley transitions. In 
this case, the Hamiltonian for the problem is given by 
(we use units such that K = 1): 



H = v F 



-kp 



where (f>^ — k x +ik y and ^ is a two component field oper- 
ator, vp is the Fermi velocity, kp is the Fermi momentum 
(related to the chemical potential fi by kp = /j,/vp), nt = 

*t + _,*£ the density operator, and Vq = ^j- the 2D 

unscreened Coulomb interaction. The non-interacting 
Hamiltonian can be easily diagonalized, giving rise to two 
different sub-bands or Dirac cones. In neutral graphene, 
the lower band is completely filled, while the upper one 
is empty. The single particle excitations in this case are 
inter-band excitations of electron-hole pairs. In the pres- 
ence of a finite chemical potential (which can be tuned us- 
ing a back gate in actual graphene devices) there are also 
intra-band electron-hole pairs in the same band. An im- 
portant point to remark here is that electrons and holes 
in different sub-bands overlap, being responsible of some 
of the most striking phenomena observed in graphene, 
like the Klein Paradox.— 

In order to describe the different approximations to 
Hamiltonian Eq. |T]), we focus on the density-density 
response function, x(<f, u), which characterizes the re- 
sponse of the system to a density probe in linear re- 
sponse theory. Its imaginary part is related, by virtue 
of the fluctuation-dissipation theorem, to the dynamic 
form factor, S$(uj) — J2 n I ( n l n ^|0) | 2 ^( w — ^no), that con- 
tains information on the relevant excitations |n) con- 
nected to the ground state |0) through the density op- 
erator (uj n o = u> n — luq is the energy of the excitation) . 

Non-interacting response function. The density re- 
sponse of free Dirac electrons is well understood in both 
undoped and doped graphene. The former was studied 
first by Gonzalez et al.£ In this paper we are primarily 
interested in its imaginary part, which gives the possible 
dissipative processes in the system (it is also connected 
to the excitations in the system): 
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graphene. In the long-wavelength region, the imaginary 
part of the polarization function is given by: 
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where Ep = vpkp is the Fermi Energy. This expression 
reflects the existence of two different particle-hole excita- 
tions with different responses: the intra-band (uj < i>Fq), 
and the inter-band (uj > 2Ep) excitations. In between 
there is a gap in the dissipative spectrum due to the fact 
that the displacement of the Fermi energy prohibits some 
inter-band transitions. Also in the doped case there is an 
infinite response at the threshold uj = vpq. 

RPA theory of interacting electrons. The RPA re- 
summation of the bubble diagram, Tl(q, uj) = — xo(<f, w) 
has been applied to a large variety of physical problems. 
As we mentioned in the introduction, when applied to 
the electron liquid it is known to capture the essential 
physics in the high density regime. Concerning undoped 
graphene, it was first employed by Gonzalez et alr^ to 
argue for a marginal Fermi liquid behavior in this sys- 
tem. In the RPA approximation the imaginary part of 
the density-density response function is: 



%XRPA(q,uj) = 



$Xo(<?,w) 



(4) 

Being proportional to the imaginary part of the bare sus- 
ceptibility, the dissipative processes allowed in the RPA 
approximation are again restricted to the region over the 
threshold. In the particle-hole continuum the real part, 
SRxoj vanishes, leading to a simplified expression for the 
RPA susceptibility: 
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We have defined the dimensionless coupling constant 
g = (g s g v /4)e 2 /vp (analogous to the fine-structure con- 
stant in quantum electrodynamics) that, in free standing 
graphene, is known to be of order g ~ 2. This approxima- 
tion is only valid for energies much smaller than the high 
energy cut-off of the theory, which is the energy band- 
width of the system A^. Actually, in the whole regime 
where the imaginary part is known to be non-zero, so is 
also the real part by the Kramers-Kronig relation: 



where g s = g v = 2 are respectively the spin and valley 
degeneracies. This function shows that free Dirac elec- 
trons have an infinite response at the threshold uj = vpq, 
and that possible dissipative processes are restricted to 
the region uj > vpq, corresponding to particle- hole con- 
tinuum. 

The case of doped graphene is slightly more compli- 
cated. It was studied by Shung in relation to interca- 
lated graphite^ and recently in Refs. [TolfTll in relation to 
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In this work, in general we will use equation ((4]) with 
non-zero real and imaginary parts. However, as we shall 
see in Section IV, including a non-zero real part is only 
relevant for energies very close to the cut-off. 

In contrast to the electron liquid, the poor screening 
properties of electrons in undoped graphene does not al- 
low to recover the Fermi Liquid picture for Dirac elec- 
trons. Even though the screening cuts off the diver- 
gence at the threshold, the Coulomb potential remains 
long ranged and there are no plasmon modes. The rele- 
vant excitations in RPA are the inter-band electron-hole 
ones, just as in the bare theory. As we shall discuss 



later, however, this picture could be uncomplete, as it 
has been proposed recently that RPA does not capture 
the whole physics in undoped graphene arguing that 
it is necessary to take into account also another sub-class 
of diagrams. 

In doped graphene this is no longer a problem, and 
the RPA approximation has been shown to capture the 
essential physics i 10 ' 11 This happens primarily due to the 
non- vanishing density of states at the Fermi level once the 
doping level is away from the Dirac Point. In the limit 
v F q << gE Fl the imaginary part of the susceptibility is: 
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UJ> 2E F 
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where luq — (9s9v£ 2 E F /2) 1 / 2 is the plasma frequency for 
graphene. In this case, the RPA shows the existence of 
a well-defined (undamped) plasmon mode in the region 
of prohibited inter-band excitations. The dispersion rela- 
tion of this plasmon is given by lu — Lu^q 1 / 2 , which has the 
same q-dependence than the 2D electron gas (as opposed 
to the density dependence of the plasma-frequency, which 
is peculiar in the case of graphene). Besides, the intra- 
band and inter-band particle-hole excitations remain in 
the dissipative spectra. 

Beyond the RPA theory in undoped graphene. As 
we have said, the validity of the RPA approximation 
in undoped graphene has been questioned.—^ The idea 
is that close to the resonance (w = v F q), there are 
non-RPA contributions which can only be neglected in 
the limit of infinite number Nf of fcrmion species (for 
graphene, we have Nf = g s g v — 4). When the condition 
Nf ~ log(v F q/\v F q — lu\) is fulfilled, it was shown that 
the correct susceptibility is given by: 



where \v is the ladder susceptibility: 



(8) 



xv(q,u) = — 



q 1 + | arcsin(a;) - (1 + ^x)^(l - x 2 ) 
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having defined x = —2= 



(9) 

ggg-lnf , v "> , ). Theimagi- 

nary part of this function is shown in Fig.[l] together with 
the other susceptibilities analyzed in the text for undoped 
graphene. The most striking feature of this susceptibil- 
ity is the new threshold of excitations that appears under 
lu = Vpq. Due to the Coulomb interaction, electron-hole 
pairs reduce their energy, giving rise to excitons. The 



new lower threshold for excitations is hence no longer 
lu = v F q, but a new one determined through the con- 
dition i = l. Also, the susceptibility shows a peak in 
this region, where the real part of the susceptibility be- 
comes positive, which can be related to a damped collec- 
tive mode. As the interaction strength is increased, the 
threshold is pulled downwards, while the collective mode 
becomes less damped. In the limit g 1 the latter is no 
longer damped, and its dispersion relation can be found 
analytically: lu q = v F q(l — e~ Nf ). In this region the con- 
tribution from excitons is negligible. Notice, however, 
that for intermediate coupling the threshold goes beyond 
the region of applicability of this susceptibility, which is 
only valid close to the resonance lu = v F q. For those 
values of the parameters a more detailed study should be 
performed. 



III. THE F-SUM RULE FOR 2D DIRAC 
ELECTRONS 

As mentioned in the introduction, exact sum rules are 
very useful as a check of different approximation schemes 
in interacting electron systems. In particular, it is known 
that in the electron liquid theory any approximation to 
the interacting Hamiltonian has to fulfill the identity: 



°° dcu q 2 N 
— w9x?,w = - — 
7r 2m 



(10) 



where N is the electron density (with volume normalized 
to one). This identity, called the f-sum rule after the 
oscillator parameters / in atomic physics, is justified for 
Schrodinger electrons, but no longer applies to the case of 
Dirac. The corresponding f-sum rule for Dirac electrons 
can be derived following the steps of Nozieres and Pines. 1 
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FIG. 1: Imaginary part of the susceptibility function for un- 
doped graphene close to the resonance, for g — 0.2. The 
thick solid line shows the result of performing the resumma- 
tion of loop and vertex diagrams. The dash-dotted line is the 
RPA approximation. The dashed line corresponds to the non 
interacting susceptibility. Notice that the former has a contri- 
bution in the region under the threshold, where Coulomb in- 
teraction gives rise to excitons and a damped collective mode 
(the highest peak) . 



The contribution from the lower filled band makes the 
sum rule cut-off dependent. This is easy to understand 
as there are inter-band particle-hole transitions for arbi- 
trary energy if we do not cut off the spectrum. More- 
over, it does not depend on the level of doping (or the 
chemical potential) since Eq. is valid also for doped 
graphene. Even more striking, the sum rule applies too 
in the case of massive Dirac fermions (i.e. in a gapped 
system) as long as the gap is much smaller than the cut- 
off: m <C Ae ■ A gap can be generated when a substrate 
breaks the graphene sublattice symmetry^ 

The f-sum rule is a statement on the conservation of 
particle number. As interactions are switched on adia- 
batically, the number of particles does not change, and 
the sum rule must be fulfilled by any susceptibility cal- 
culated approximately from the interacting Hamiltonian. 
The effect of interactions is usually a rearrangement of 
the contribution to the f-sum rule by the different ex- 
citations. This is immediately related to the spectral 
weight (|("|wt|0)| 2 ) ea; of the excitation, a quantity that 
is relevant in dissipative processes^ and other physical 
properties such as van der Waals forces between neutral 
systems J£ 



Details of the derivation are given in the App. [A] The 
result is: 



gs9vq 2 A E 
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where Ae is the energy cut-off (i.e. of the order of the 
band- width). The f-sum rule was studied earlier in the 
context of the true 3D relativistic Dirac problem^ with 
infinite band- width (no lattice). In that case, the sum 
rule gives an infinite contribution that is considered part 
of the "vacuum" energy and hence unmeasurable. In 
a system with a finite band-width, which is the case 
of graphene, such a cut-off dependent contribution rep- 
resents a real response of the system (the inter-band 
electron- hole excitations), and hence is non- negligible. 
This issue was raised originally in Ref.— , where a rel- 
ativistic generalization of the f-sum rule is discussed too. 
Notice that, as shown in the appendix, in order to de- 
duce a non-zero f-sum rule for Dirac electrons we must 
proceed carefully, since conventional general definitions 
of the f-sum rule do not apply immediately due to the 
peculiar structure of the Dirac Hamiltonian. 

The f-sum rule for the imaginary part of the suscepti- 
bility can be related to a f-sum rule for the optical con- 
ductivity straight-forwardly. In a recent worki£, the lat- 
ter has been studied for the whole band of graphene, by 
using a tight-binding description. Though this approach 
is more realistic when applied to graphene, the f-sum 
rule for the continuous case provides more insights when 
studying models of interacting electrons in the low-energy 
regime, which is our interest in the present work. 

The f-sum rule for Dirac electrons have many pecu- 
liarities that are not present in the usual Fermi liquid. 



IV. TRANSFER OF SPECTRAL WEIGHT 

We now turn to the application of the exact f-sum rule, 
derived in the last section, to particular approximations 
for the Coulomb-interacting Dirac problem, stressing the 
(re)distribution of the spectral weight among different ex- 
citations. Firstly, we study the case of free electrons, in 
both doped and undoped graphene. Then we discuss the 
more controversial issue of graphene with interactions. 
Here, the discussion of Section II on the different ap- 
proximations to the problem will be especially relevant. 
From here on we concentrate on the particular case of 
Nf = 4 fermion species, which is the case of graphene. 



A. Free Dirac electrons 

If graphene is undoped, i.e. the Fermi level is located 
at the Dirac Point, the susceptibility is given by Eq. ^j. 
This expression saturates the f-sum rule for Dirac elec- 
trons, as expected, meaning that the contribution from 
inter-band particle-hole excitations has a q 2 dependence 
and that the number of those excitations is only limited 
by the high energy cut-off. 

When the level of doping is changed away from half fill- 
ing, we have two different particle- hole excitations: intra- 
band and inter-band. As the f-sum rule does not depend 
on doping, the spectral weight must be distributed among 
those. By using Eq. ([3]), in the long- wavelength limit, we 
find a contribution from inter-band transitions: 
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On the other hand, the contribution from the intra-band 
excitations is: 



2k f 
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q 2 E F 



(13) 



meaning that there is an exact transfer of spectral weight 
between the two types of excitations that is proportional 
to the level of doping. When the Fermi level is changed, 
some inter-band transitions are prohibited by Pauli ex- 
clusion, and they no longer saturate the f-sum rule. Both 
excitations have a contribution proportional to q 2 . How- 
ever, the contribution from intra-band excitations is not 
cut-off dependent, but density (or doping) dependent. As 
we shall see when connecting the interactions, those ex- 
citations resemble those of the conventional Fermi liquid, 
where the contribution to the f-sum rule due to particle- 
hole excitations does depend on the density of electrons. 



inter-band excitations is again given by Eq. (fT2|) . which 
states that these excitations have lost a contribution to 
the f-sum rule given by E F q 2 /2. As opposed to the non- 
interacting theory, this contribution is not transferred to 
the intra-band particle-hole transitions. In this case their 
contribution is: 
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which is proportional to g 4 and thus much smaller in the 
long-wavelength limit than the loss from the inter-band 
transitions. As in the case of the Fermi liquid, it is the 
plasmon mode which absorbs most of the spectral weight: 
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B. Interacting Dirac electrons 

RPA in undoped graphene. As we have seen in Section 
II, the simplest approximation we can employ to study 
the role of interactions is the RPA. For an undoped sheet, 
we use equations (@| and © and find that the RPA 
approximation fulfills the f-sum rule. The main issue 
related to the RPA comes when analyzing the spectral 
weight. In principle, in undoped graphene at the RPA 
level, there are no new low-energy excitations besides the 
particle-hole. However, it can be seen that the low-energy 
sector loses some spectral weight. If we take an interme- 
diate energy scale, say A/, such that v F q « A; « A F 
and integrate over energy in this range, we get: 



q 2 Ai 



—gv F q 3 . (14) 



Thus there is a cut-off independent loss of spectral weight 
in the low-energy sector. This is an unconventional result 
in the sense that the usual RPA approximation for Fermi 
liquids only rearranges the low-energy spectral weight, 
and does not couple to the high energy sector of the the- 
ory. On the contrary, in undoped graphene it seems to 
be a spectral transfer to the high energy sector. In par- 
ticular, the spectral weight is distributed close to the 
cut-off energy scale, where we find a new resonance and 
even a plasmon condition that is fulfilled. Of course, the 
particular kind of high energy excitations cannot be de- 
scribed by our effective theory, which only applies to the 
low-energy regime. The transfer is proportional to the 
interaction strength, <?, and it depends on the momen- 
tum as q 3 , instead of the typical q 2 dependence of single 
particle excitations, being negligible at leading order in 
the long- wavelength limit. 

RPA in doped graphene. In the doped case, the RPA 
susceptibility is given by Eq. flT]), which predicts three 
relevant excitations: inter-band and intra-band particle- 
hole, and a plasmon mode. The contribution from the 



Hence, in the asymptotic limit of q —> 0, the spectral 
weight of the intra-band excitations is exactly absorbed 
by the plasmon. 

In the interacting system the response is dominated by 
the plasmon, instead of the intra-band particle - hole ex- 
citations. Notice, however, that graphene, due to its par- 
ticular two-band structure, is different from the electron 
gas in the sense that the main contribution comes always 
from the inter-band transitions, which essentially satu- 
rates the f-sum rule. Despite this fact, as will be shown 
in next section, the existence of this collective mode can 
have measurable consequences. 

Far above the energy scale E F , we recover the same 
polarizability that for undoped graphene, as expected 
as now the only relevant excitations are the intra-band. 
Hence, a transfer of spectral weight to the high energy 
sector is also observed. Notice that, at the same order 
of this transfer, q 3 , further rearrangements of spectral 
weight occur in the low-energy sector, as can be shown 
by computing next-order corrections to expression ([7]). 

Beyond RPA in undoped graphene. As mentioned 
above, in doped graphene the RPA approximation 
seems to work well, so we concentrate now in undoped 
graphene. As we have seen, when vertex corrections are 
taken into account, the polarization function close to the 
threshold w — v F q is given by Eqs. §8§ and ©■ As 
this approximation does not describe the whole energy 
range, we cannot verify immediately the validity of the 
f-sum rule. However, we know that for high energies the 
RPA approximation should be valid, so we expect a loss 
of spectral weight similar to the one found before. The 
most relevant question is how the spectral weight is dis- 
tributed among the new excitations that arise. 

We would like to study the problem of spectral weight 
transfer when vertex corrections are added. As its con- 
tribution to the f-sum rule is related to the role they play 
in the response of the system, we could expect that new 
excitations can give deviations from the free electrons pic- 
ture. However, the new modes found in Ref. HI can only 
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FIG. 2: Contributions to the f-sum rule in the low-energy 
region for undoped graphene. The integrals are done with a 
high energy cut-off such that ui < l.lvpq. The dotted line 
corresponds to the non-interacting theory. The dash-dotted 
line shows the RPA contribution. The solid line is the non- 
RPA contribution, sum of the excitonic contribution (under 
the threshold, u> < vf<i) and the particle-hole contribution 
(over the threshold, u> > VFq)- 

be calculated reliably close to u> = vpq, and one cannot 
integrate the susceptibility given in Eq. ^ in the entire 
energy range. Nevertheless, we can concentrate on the 
transfer of spectral weight close to lu = vfQ by defining 
the following quantity: 

4 f Aq 

S P (g,C) = 3/ duj uj ^x{Q, w ) . (17) 

vf<t Jo 

where A q — CvFq and C is a constant that has to be 
chosen appropriately. If we define s = lj/vfQ we find: 

S P (g, C) = [ C d SS (4v F ^x(q, s)/q) (18) 
Jo 

which is momentum independent but cut-off and interac- 
tion strength dependent. We find that the approximation 
of Eq. works well for C — 1 ~ 1CP 1 (in such a way 
that Nf = 4 ~ logoff?/ (A« — v Fq)j) and this is the value 
we use from now on. Readily, we notice that all the con- 
tributions to the f-sum rule that come from this region 
are of order q 3 , and hence all the spectral rearrangement 
because of interactions will be of subleading order. 

We analyze in the following the dependence on the 
interaction strength g. In Fig. [2] we plot the various con- 
tributions to S p (g, 1.1). This Figure shows a similar be- 
havior for the electron-hole background in both RPA and 
non-RPA approximations: as the interaction strength is 
increased, the spectral weight covered in the integration 
is smaller. Notice that this should not be related only 
to the transfer of spectral weight to the high-energy sec- 
tor (which is linear in g), but mostly to the fact that 
the maximum of the imaginary part of the susceptibil- 
ity is shifted away from the threshold as the interaction 



grows. More interesting is the behavior of the excitonic 
part: it shows a maximum for g ~ 0.2, and then starts 
to decrease. 

In order to understand this behavior, we notice that 
the contribution to the f-sum rule from the plasmon can 
be easily worked out in the limit g 1, where the plas- 
mon gets well defined: 

/ d WW (9f X ) pI = -^[^(l-e-^)] J (19) 
Jpi 4 9 

and where we have used the analytical plasmon disper- 
sion relation for generic number of fermion species (re- 
member that in the case of graphene we have Nf = 
4), ui(q) — vpq(l ~ e~ Nf ), and its spectral weight, 
|(n|^|0)|pj = 2vFq 2 /ge~ Nf . We see that the plasmon 
contribution to the f-sum rule decreases as the interac- 
tion strength grows. On the other hand, although the 
excitonic domain gets broader (in thus out of the region 
of validity of our approximation) , its contribution is also 
negligible for large interactions. Hence, the appearance 
of a maximum in the contribution to the f-sum rule under 
the threshold can be understood as an interplay between 
excitonic response and plasmon response. 

Anyway, our results show that there is no important 
transfer of spectral weight from the particle - hole ex- 
citations to the new modes predicted in the context of 
vertex corrections, as they are of order q 3 . As in prin- 
ciple this approximation only rearranges spectral weight 
close to the threshold, we expect a flow of spectral weight 
to the high energy sector of the theory as the interaction 
strength is increased, as occurred in the RPA approxi- 
mation. 



V. ANALYSIS OF THE RESULTS 

The different contributions to the f-sum rule have been 
summarized in Table |TJ The Dirac liquid shows partic- 
ular features from the point of view of spectral weight 
transfer and contributions to the f-sum rule. The results 
provide many insights into the role played by Coulomb 
interactions in graphene. 

The f-sum rule is essentially saturated by the particle- 
hole excitations in the long-wavelenght limit. This is true 
for small doping levels compared to the cut-off, as the 
intra-band excitations have the largest contribution. In 
the case of doped graphene, Coulomb interactions give 
rise to a collective mode which absorbs to leading order 
all the spectral weight of inter-band excitations. How- 
ever, this does not remain true in the case of undoped 
graphene, where all the spectral (re)arrangement due to 
the Coulomb interaction is always at subleading order, 
q 3 . New excitations which could be arising close to the 
threshold contribute much less to the f-sum rule than the 
particle - hole ones, even when the interactions strength 
is infinitely large. This effect must be understood not 
only because of the subleading dependence of those exci- 
tations, but because of the spectral weight transfer to the 
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Free 


Free 


RPA 


RPA 


Beyond RPA 




Undoped 


Doped 


Undoped 


Doped 


Undoped 


Inter-band 


A E q 2 


(As - 2E F )q 2 


A E q 2 


(A E - 2E F )q 2 


A E q 2 


Intra-band 




2E F q 2 




vW/g 2 E F 




Plasmon /excitons 








2E F q 2 


v F q 3 /g 


High energy sector 






gv F q 3 


gv F q 3 


gv F q 3 



TABLE I: Contributions to the f-sum rule from the different 
excitations present in the theories of non-interacting and in- 
teracting electrons in graphene studied in this work, both for 
doped and undoped graphene. 



high energy sector we have found in every approximation 
studied for the interacting Dirac liquid. Hence, the latter 
should be a remarkable feature of this system. 

An apparently related behavior has been described in 
high-Tc literature, where spectral weight transfer be- 
tween different energy scales has been reported in several 
works . 19 i 20 In the context of strongly correlated systems, 
a spectral weight transfer from the high energy sector 
to the low-energy one has been related to possible issues 
arising in the definition of the low-energy theory— In 
this sense, studies concerning the whole graphene band- 
structure could help to clarify this issued 

When turning to graphene, the Dirac liquid is only 
a low-energy approximation to the electronic structure, 
and the f-sum rule derived here does not cover the whole 
band. However, from our results we expect in this case 
a transfer of spectral weight beyond the artificial cut-off 
that we introduced to delimitate the continuous descrip- 
tion. Actually, some related effect has been observed in 
the Coulomb impurity problem in graphene^ where a 
bound state appears beyond the band, i.e., in the high 
energy sector. 

The relative lack of importance of Coulomb effects in 
graphene is in agreement with the experimental obser- 
vation that electronic carriers are very well described 
by free Dirac fcrmions. This is specially true for un- 
doped graphene, but also applies to doped graphene, 
where dissipative processes would be dominated by the 
excitation of the incoherent particle - hole background. 
Systematic experimental studies, however, should give 



a clear trace of the plasmon mode in doped graphene, 
which actually seems to have been observed in ARPES 
experiment s 23 ! 24 ' 25 and could be in principle detected in 
inelastic x-ray scattering experiments, where information 
about the dynamical structure factor can be extracted. 
Another physical quantity where plasmons could be im- 
portant is the van der Waals force between graphene and 
other neutral systems, i.e., another graphene layer or a 
substrate^. In perturbation theory, the first correction 
to the energy of the system can be shown to be: 

E m =-J2 [ duV q 2 (z) XG (q, iuj) X G(q, iw) (20) 
q J 

In undoped graphene, despite the new excitations pre- 
dicted, this correction can be shown to be = 

— Vi(^-) 3 in leading order^ where Zq and Vj are, respec- 
tively, typical distance and energy scales of the problem, 
and z is the distance between the two graphene sheets. 
On the contrary, in doped graphene, the response of the 
plasmon gives a new contribution: 

E^ = -V^f-V pl {^f' 2 (21) 

For an intermediate regime of distances, the dominant 
interaction is the one coming from the response of particle 

- hole excitations, as V* 3> V v \. The leading contribution, 
however, should be the one coming from the plasmon, as 
it is the one which decays more slowly. 



VI. CONCLUSIONS 



We have derived the f-sum rule for Dirac electrons and 
applied it to study the spectral transfer in different ap- 
proximations for interacting electrons in graphene. The 
f-sum rule, being an statement of the particle number 
conservation, is not interaction dependent and hence has 
to be fulfilled by any theory of the interacting Dirac liq- 
uid. Hence, it provides an insightful tool to study the 
effects of electron-electron interactions in the low-energy 



regime of graphene samples. 

We have shown that the RPA theory for undoped and 
doped graphene fulfills this identity, even though we find 
that its behavior turns out to be quite unconventional. 
Instead of having a spectral weight transfer among states 
around the same energy scale, the transfer involves two 
different scales, with one being on the order of the cut-off. 
When vertex corrections are included, this behavior is 
not changed, being a feature of the Dirac liquid. Despite 
the f-sum rule derived here only applies to low-energy 
graphene, we expect a similar transfer of spectral weight 
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to occur once the whole band is taken into account. 

Besides, we have studied the relative importance of 
the different excitations predicted close to the Dirac 
point in graphene. The f-sum rule is essentially satu- 
rated by the inter-band particle- hole excitations, though 
in doped graphene the collective plasmon acquires some 
importance which could have measurable consequences in 
ARPES or x-ray inelastic scattering experiments. This 
is not the case in undoped graphene, where any new ex- 
citation coming from Coulomb interaction is essentially 
negligible compared to the electron-hole ones, as far as 
the spectral weight is concerned. Surprisingly, this fea- 
ture remains true even for the collective mode predicted 
when vertex corrections are taken into account. Hence, 
its possible experimental observation, as well as the ob- 
servation of other excitations different from the electron- 
hole ones in undoped graphene, can be technically chal- 
lenging. As far as the measurements imply a certain av- 
erage of the spectral weight, we expect that undoped 
graphene will respond to experimental probes essentially 
as a non-interacting system. 
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APPENDIX A: DERIVATION OF THE F-SUM 
RULE 

The f-sum rule can be derived by relating the integral 
over the first momentum of the structure function and 
the following commutator>i 

/>oo 

(0|[[n f ,tf],4]|0)=2 / dwuSfc). (Al) 
Jo 

For what concerns us, it is more useful to express this 
identity in terms of the density-density correlation func- 
tion, x(<Z»: 

2 f°° 

(0\[[nf,H],nt]\0) = — / dwwQft^w). (A2) 
q n Jo 



1. Undoped graphene 

In graphene, at low energies, if we focus on a single 
valley and neglect the spin degree of freedom, the Hamil- 



tonian is: 

*°=-E*:-( fe+i ^-*s)% (A3) 

where \l/t = (ag, bj^) is a spinor, being a and b operators 
referring to the two different sublattices in the unit cell 
of the honeycomb lattice. In the same basis, the density 
operator is given by: 

nt=y*t % = V(at ar + bl 6A (A4) 

1 k+q * k+q * k+q k> \ ' 

k k 

The Coulomb Interaction has the usual form Hi nt = 
^qVq n q n q ( U P t° a term that is proportional to the num- 
ber of particles in the system that is not relevant for this 
discussion), so it is simple to show that the condition 
[n q ,Hint] = is realized. So wee see clearly the advan- 
tage of the f-sum rule: it can be calculated within the 
non-interacting theory, but it must be satisfied also by 
the interacting one. 

With these expressions, the commutators in Eq. (|A2|) 
can be evaluated. We notice that the first commutator 
is nothing but the particle conservation equation written 
in the momentum space: 

[n^H] =q- f ? , (A5) 

where J$ = vf Y]g ^to^g^- is the velocity (current) op- 
erator. So the double commutator reads: 

% 

Now we should proceed carefully with this result. In prin- 
ciple, our free theory describes massless electrons with 
an unbounded linear dispersion relation. But when han- 
dling operators that are defined in an unbounded region, 
we are not allowed, for instance, to simply state that 
^ fc G(k + q) = ^2 k G(k). The same issue can be found 
in the theory of the one dimensional electron liquid when 
calculating the commutator [n q , n q \ , which turns out to 
be nonzero, giving rise to the anomalous commutator 
problem. 2 Once we work with unbounded operators, we 
need to refer the calculation to bounded quantities, which 
usually are defined with respect to the ground state value. 
Then, we define normal ordered operators: 

: G(k) : = G(k) - (0|G(fc)|0). (A7) 

As now the normal ordered operator is bounded, it sat- 
isfies J2k ■ + 1) ■ ~ J2k ■ ^ W : ■ Applying this rule to 
our commutator, we get: 

4 = 

k 

(A8) 
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where the ground state for undoped graphene consists on 
the lower band completely filled (Dirac sea) and the up- 
per band completely empty. The result involves the dif- 
ference between two infinite sums, something not defined 
a priori. To compute it, we need to regularize somehow 
the sums, for instance, by using an ultraviolet cut-off. 
Besides, we need to switch to the diagonal basis. As 
usual, this is done with the unitary transformation that 
diagonalizes Hamiltonian (IA3[) : 

U^^( er f e ~t), (A9) 
which leads to the following relation: 

\-i sm(6$ - 9 Uq -, ) - cos(0 9 - - Uq -, ))' 

Once we have the cut-off we can shift the sums: 

(0\[[n ? ,H},n^\0) = - Wf [^cos(%) - £ cos (%)" > 

kei ken 

(All) 

where only the contribution from the lower band survives, 
and the regions of summation are shown in Fig. [3l For 



2. Doped graphene 

The effect of doping the graphene sheet translates in a 
nonzero chemical potential. The ground state no longer 
has a particle-hole symmetry. Therefore, the modifica- 
tions to the f-sum rule calculation of the last section are: 
(1) The Hamiltonian must be replaced by H — fiN, with 
fi being the chemical potential. (2) The ground state has 
a contribution from electrons of the upper band, for pos- 
itive chemical potential, or from holes in the lower band, 
for negative chemical potential. 

First we discuss (1). We have N = n^, which leads to 
[n q , N] = 0, thus giving no new contribution. Therefore, 
at the level of commutators the result is again Eq. (jlip . 
From (2) we could expect a contribution from the elec- 
tron (holes) in the ground state, above (below) the Dirac 
point. However, this contribution is no longer unbounded 
(the operators are only nonzero under (over) the Fermi 
momentum), and now we can shift the operators: 

^cos(%-^ +? )(0|ct +fC£+f |0) 

- cos(%- £ )(O|ct C£ |O) = 0. (A13) 

Therefore, there is no contribution coming from the 
ground state evaluation of the commutators, and the re- 
sult is the same than the undoped one, Eq. (fTTj) . 




FIG. 

rule 



3: Regions subtracted in the calculation of the f-sum 



a large momentum space cut-off A = Ae/(2vf), the cal- 
culation of the difference can be carried out easily. Both 
regions give the same contribution but with opposite sign 
due to the cosine term. By going to the continuum limit 
of the sum we find: 



q 2 A E 



(A12) 



Reinstating the degeneracy of spin and valley gives rise to 
an extra multiplicative factor g v g s — 4. The final result 
for the f-sum rule is Eq. (fTTj) : 



gsg v q 2 k E 
16 



3. Massive Dirac electrons 

The case of massive electrons requires more attention. 
The Dirac Hamiltonian in this case reads: 



m 

k * vp(k x + ik 



v F (k x 




*£, (A14) 



which is the same as Eq. (| A3|) with the additional term 
mJ2k ^k a ^k- Again, modifications to the f-sum rule 
could arise from the new term in the commutator or 
from the final ground state evaluation. The first con- 
tribution can be readily seen to be zero: [n q ,H m ] = 0, 
due to the cancellation of the sublattice contributions in- 
dependently. The second modification is more subtle, as 
requires the diagonalization of Hamiltonian in Eq. (|A14[) . 
The result gives an hyperbolic dispersion relation of the 
form Ek = \/vpk 2 + m 2 . Applying it to the commuta- 
tors which give the f-sum rule (following the same lines 
than above) we get: 



n ? ,fl],nj|0) = 



-V F 



rr^ y/ (vFk) 2 + m 2 
vpk 



cos(%) 



E 

ken 



\J (v F k) 2 + m 



2 COf H^ 



(A15) 
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The result again involves the subtraction of two regions, 
which is nonzero only for momenta close to the cut-off, 
as we showed in Fig. [3] If m << we can expand the 
prefactor to leading order in k, recovering the massless 
result. This is to be expected since the mass term is irrel- 
evant in this region. The result for the massive case turns 



out to be the same than the massless one of Eq. (fTTj) , so 
the existence of a gap does not change the nature of the 
sum rule. At this point we remark that the particular 
form of the f-sum rule for Dirac electrons comes mainly 
from the asymptotic linear spectrum and the existence 
of an unbounded spectrum. 
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